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Abstract 

Is there a functional role for quantum mechanics or coherent quantum effects in biological processes? While this 
question is as old as quantum theory, only recently have measurements on biological systems on ultra-fast time-scales 
shed light on a possible answer. In this review we give an overview of the two main candidates for biological systems 
which may harness such functional quantum effects: photosynthesis and magnetoreception. We discuss some of 
the latest evidence both for and against room temperature quantum coherence, and consider whether there is truly 
a functional role for coherence in these biological mechanisms. Finally, we give a brief overview of some more 
speculative examples of functional quantum biology including the sense of smell, long-range quantum tunneling in 
proteins, biological photoreceptors, and the flow of ions across a cell membrane. 
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1. Introduction 

The role of quantum mechanics in biological processes has a long history 1110, Eh- After all, biological molecules 
are constructed of atoms, and bound together by forces of a 'quantum origin' . However the role of quantum coher- 
ence (and more recently entanglement) on time or energy scales which play a direct role in biological function has 
remained controversial. Recent experiments |j4l |5l la tZl ISl on photosynthetic 'Light Harvesting Complexes' (LHC) 
and their constituents (e.g., the Fenna-Matthews-Olson (FMO) pigment-protein complex in green sulfur bacteria) 
have suggested that quantum coherence may play a role in one of the most fundamental and important of biological 
processes: energy transport and energy conversion. 

LHC complexes are arrangements of pigments (most importantly chlorophylls) and protein molecules which func- 
tion as a light gathering 'antenna' to absorb photons and become electronically excited 0[l3l. This excitation is then 
passed to a reaction center, where it is converted into useful chemical energy. One of the most interesting features 
of this biological machine is the highly efficient transfer mechanism which takes the electronic excitation through the 
LHC to the reaction center with almost unity quantum yield. In green sulfur bacteria, part of this chain connecting the 
antenna (which in green bacteria is a large chlorosome) to the reaction center is composed of the Fenna-Matthews- 
Olson (FMO) pigment-protein complexes ITll; see Fig. 1 for an explanatory schematic. It is this complex which has 
been most widely studied, and will form part of the focus of this review. 

In photosynthetic light harvesting quantum coherence can have different meanings in different circumstances. 
Quantum coherence can refer to quantum superpositions of localized molecular excitations that occur naturally be- 
cause electronic couplings between molecular excitations lead to delocalized eigenstates (a.k.a. excitons). Sometimes, 
coherence in light harvesting alludes to the coherent wave-like dynamics of energy transfer, which actually reflects 
the superposition of excitonic eigenstates. In the former case, the coherence is represented in the molecular site basis, 
whereas in the latter case, the coherence is represented in the delocalized exciton basis. Note that the site basis and the 
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exciton basis are special because they are related to the spatial arrangement of chromophores and energy eigenstates 
of the Hamiltonian, respectively. Although both types of coherence effects play important roles in photosynthetic 
light harvesting, they must be discussed separately. Quantum coherence manifested in the delocalized eigenstates 
of photoexcitations in photosynthetic complexes plays a fundamental role in spectral properties, energy tuning, and 
energy transfer dynamics of photosynthetic light harvesting 1 12, TslIiSl- The effects of excitonic coherence are more 
difficult to analyze and have been the subject of intensive research in the past few years. Spectroscopic measure- 
ments (termed two-dimensional electronic spectroscopy) on the FMO complex at liquid-nitrogen (77K) 141] and room 
temperatures have shown time-dependent oscillations, presumed to be quantum beating, in the amplitudes of the 
spectral signals, which matches the predictions of quantum theory. Inspired by these observations, it has been pro- 
posed that a intricate interplay of quantum coherent excitation transfer and environmental dephasing help increase the 
efficiency of the energy transport process 

Apart from the example of photosynthesis it is also believed that other functional biological mechanisms may 
display, or rely upon, quantum effects. For example, the second part of this review will focus on the singlet and triplet 
states of spatially-separated electron spins in radical pairs (donor-acceptor molecules), whichh are hypothesized to 
play a role in avian magnetoreception |2^ 21, 22|. If true, this would represent a functional piece of biological 
"quantum hardware" that could not function in a classical worl d II23I1. More speculative ex amp les we will describe 
more briefly include long-range quantum tunneling in proteins |24l 12511 . the sense of smell 12611 . biological photore- 
ceptors 1I27II28II . and the flow of ions across a cell membrane 112 911 . 



1.1. Functionality defined 

Our goal here is to give a brief overview of these different biological mechanisms, summarize some of the theo- 
retical models, and highlight some of the experimental and theoretical evidence both for and against a functional role 
for quantum effects in biological systems. By "functional" we imply a role where the presence of coherent quantum 
dynamics achieves something either more efficiently, or otherwise impossible, than could be achieved by a classical 
mechanism alone. This concept of nature taking advantage of quantum mechanics is an inspiring notion, but the ev- 
idence for and against it must be examined carefully. We thus discuss only experimentally-verifiable systems in this 
review. This topic is growing at a phenomenal pace, and we can only summarize the main points for each system we 
consider 



7.7.7. Photosynthesis 

In the example of energy transfer in photosynthetic complexes, most of the quantum effects observed so far 
have been at liquid-Nitrogen temperatures (77 K) though recent evidence has arisen of persistence quantum 
coherence at room temperature, both in FMO JtI and in other components of certain light-harvesting complexes 
10]. To compliment this a variety of theoretical techniques have been used to show that a quantum model of the 
FMO complex implies a high transport efficiency even at room temperature llisl Il6i Il7l boi bill . As mentioned. 



this amalgamation of experimental observation and theoretical modeling suggests that a combination of quantum 
effects and environmental 'stimulation' leads to a high efficiency in the rate of energy transport in the FMO complex. 
However, alternative points of view still exist in the literature ll32ll . and more work remains to be done to analyze to 
what extent this environment-assisted transport is really 'quantum'. For example, in principle one must eliminate all 
possible classical models ll33ll and purely thermal/environmental effects 134] that could produce similar efficiencies 
before unambiguously stating that high efficiency transport itself is an example of functional quantum biology. In 
addition, it is not yet clear if the dynamics observed in-vitro (in experiment) also occur in-vivo (in nature), due 
to uncertainty about the nature of the energy transport from the antenna (chlorosome and baseplate) to the FMO 
complex. We will attempt to summarize all of these issues in this article. 



7.2. Avian magnetoreception 

In contrast to photosynthesis, the functionality of the proposed radical pair model for avian magnetoreception is 
entirely dependent on quantum mechanics 12311 . In this sense it is different from the possibly more ambiguous role 
of quantum mechanics in photosynthesis. For example, the proposed mechanism of radical-pair magnetoreception 
cannot function at all without a large degree of spin-spin entanglement and coherence 12311 : if environmental noise 
is too strong (an effective classical limit of the model) then birds could not navigate via this proposed mechanism. 
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In addition, while ubiquitous in molecular and atomic physics, there is essentially also no classical analogue to the 
singlet/triplet states of two coupled spins, which is the fundamental element of this model which acts as a sensor of 
the Earth's magnetic field. 

The open question in this case then is whether such a proposed radical pair really exists and functions as pre- 
dicted. There is some suggestion it may reside in cryptochromes (a light-sensitive cell found in the eye). However, 
experiments on a range of possible radical-pair molecules in laboratory conditions have not yet found a radical-pair 
molecule with the desired properties. Furthermore, evidence must also be found for whether birds really have the 
biological circuitry in place to process the signal that radical-pair magnetoreception would generate. Later we will 
summarize the details of this mechanism, explain why it is a strong contender for functional quantum biology, and 
outline the evidence supporting it as a mechanism for magnetoreception. 

2. Photosynthesis 



Biological system Results 
Photosynthesis 

Cyrogenic temperature quantum coherence 
Ambient/room temperature quantum coherence (FMO) 
Ambient/room temperature quantum coherence (Algae) 
Environment assisted transport 
Entanglement, Leggett-Garg 
Alternative views 

Radical Pair Magnetoreception 

Early proposals and evidence 
Mathematical models 
Indirect evidence (light dependance, magnetic field) 
Experiments on Radical pairs 
Monarch butterflies, Elys, Humans 

Other examples 

Olfaction [26,54] 
Ion channels [291 
Vision [27,55J 
Long-range electron transfer [24. ^5J 

Enzyme catalysis [56 , 571 

Table 1 : Summary of a selection of the main experimental and theoretical works on functional quantum biology. 

Photosynthesis is one of the most important photochemical processes on Earth. As outlined in the introduction, 
the primary photosynthetic apparatus (which is typically called a photosynthetic unit, or PSU) can be roughly divided 
into several parts: a light-harvesting complex antenna, and a reaction center (RC), as shown in Fig. 1 for the example 
of green sulfur bacteria. The function of the FMO complex in green sulfur bacteria is to transport excitations from 
the sunlight-harvesting LHC antenna (chlorosome) to the reaction center (RC), where it initiates a charge-separation 
process that generates chemical potential for biochemistry. 

The FMO complex consists of eight bacteriochlorophyll-a (BChl-a) molecules which are bound to a surrounding 
protein scaffolding. Note that the eighth BChl was not discovered in structural models of the FMO complexes until 
recently, therefore most studies on the FMO complex so far only considered a model with seven BChls. As we will 
describe later, in a simplified picture, each BChl molecule can be excited from its ground state into its first singlet 
excited state, forming a molecular exciton. Electronic couplings between molecular excitons then enable excitation 
energy transfer between BChl molecules. 

After a photon is caught by the chlorosome antenna, it is transferred through the baseplate to the FMO complex 
(see Fig. 1). This excitation is then passed from one BChl molecule to the next, until it reaches the molecule, or site, 
closest to the reaction center It then irreversibly enters the reaction center, and ignites the charge-separation process. 
In general terms, this apparatus functions in a similar way across a broad range of photosynthetic organisms, though 



15, 4J 

U5,i6,3S,i9J 

[36,37] 
138,331 



□£], [4r] 
[42,43,44,45,46,47] 
[48, 49, 50, ST] 
[52,531 
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the size and configuration of the antenna and component systems can differ greatly 115 811 (e.g., see Hu et al 115 911 for a 
review of the PSU in purple bacteria). 

The most well-known way to describe this excitation-transfer process is the Forster model lIHol . In this model 
the dipole-transition coupling between two chromophores (sites) gives rise to rates for excitation (exciton) hopping 
between interacting sites. The Forster model is a successful one, but it neglects quantum coherences between different 
sites. Experiments have shown that the exciton can move coherently among several chromophores, leading to elec- 
tronic quantum beatin g 116 ill, a nd perhaps also coherent collective phenomena like the well-known optical example of 



superradiance 1621 l63l l64 l65l 



In addition, some other components of LHCs in other species suggest that coherence and long-range entanglement 
may exist in larger structures. For example, the LHC of certain types of purple bacteria contain two types of ring- 
like antenna, LHl and LH2 [s^j, which are substantially larger than the FMO complex. Spectroscopic studies on 
these systems have suggested that photoexcitations may delocalize among 4-5 BChl sites in the strongly coupled 



rings ll62ll63Ll65 1. 

Another interesting issue is that the energy transport through the FMO complex is very fast (lOO's of femto- 
seconds), and the efficiency of many LHC systems is very high. The quantum yield of light-to-charge conversion in 
some photosynthetic units ifioll can be up to 95%, and explaining how the excitations navigate the energy landscape 
of the LHC so successfully is non-trivial. In addition, the need to navigate this landscape quickly must be understood 
in terms of losses that occur during the transport through the LHC system and its components, and not in terms of the 
overall energy processing rate of the entire PSU (LHC and RC systems), which is relatively slow in comparison. For 



example, the 'reset time' or 'turnover' rate of the RC in purple bacteria 115911 is of the order of 1 kHz, and the energy 
absorption rate of a single BChl in an antenna in bright conditions is 10 Hz. It is thought that the large number of 
BChl in the antenna in totality provide enough energy to optimally use this 1 kHz turnover rate of the RC. Thus there 
is a large separation of time-scales. The fast transport through the FMO complex, for example, is needed to "beat" 
the rate at which excitations are lost due to fluorescence relaxation, not because of the overall "clock rate" of the light 
harvesting complex. Also, one should note that the "energy" efficiency of photosynthetic bio-mass as, e.g., a fuel 
source, is relatively low even compared with photovoltaic efficiency 166116711 . because the down-stream biochemical 
reactions that turn chemical potentials into biomass have extremely low energy efficiency. 

To gain a deeper understanding of these quantum phenomena, powerful spectroscopic tools have been developed, 
which provide opportunities to study the dynamics of the excitations in these pigment-protein complexes at the time 
scales of femto-seconds Il65u68il . For example, two-dimensional electronic spectroscopy is a powerful technique that 
probes the couplings and the dynamics of energy flow on a two-dimensional (2D) map in the frequency domain that 
allows direct observation of "coherence" between electronic excitations. A full understanding of exactly how much 
information can be extracted with such techniques is still being discussed and researched 1i69l 70[ 71 j. 



2.1. Experimental signatures of quantumness 

As mentioned, the advent of powerful experimental tools have opened-up the ability to probe quantum effects in 
photosynthetic systems with unprecedented sensitivity. In particular, 2D electronic spectroscopy provides a unique 
tool that is specifically sensitive to electronic coherence in excitonic systems |74, 3 25, 2^. A 2D experiment is 
a four-wave mixing process, in which three laser pulses interact with the sample to create an electronic polarization 
that generates the signal (Fig. |2^). In the experiment, two periods of time delay can be controlled in the apparatus; 
the time delay between the first and second pulses (r, coherence time) and time delay between the second and third 
pulses (r, population time). At fixed r and T, the signal field emitted in the photon-echo phase-matching direction 
ks - -k\ + k2 + kj, is combined with an attenuated local-oscillator pulse for heterodyne-detection and frequency 
resolved to obtain a spectrum, which can be regarded as the Fourier transform of the oscillating signal fields with 
respect to the time delay between the third pulse and the signal (f, re-phasing time). The measurement can be repeated 
at varied t and Fourier transformed with respect to t to obtain a 2D spectrum in frequency domain {lOj and o),) at a 



fixed population time T 11751 17411 



Because 2D electronic spectroscopy records the signal at the level of the field rather than the intensity, it is 
sensitive to the quantum phase evolution of the electronic system during the population time. A broadband pulse can 
then interact with multiple exciton states to produce superpositions (coherences) of them, and the induced coherence 
in the exciton basis then undergoes oscillatory phase evolution that leads to beating signals as a function of time 
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reaction center 



Figure 1: (Color online) (a) Schematic diagram (from f72ll ) of roughly 200,000 BChl-c molecules (green rhombuses) encased in a proteinlipid 
structure forni the chlorosome antenna. This is much larger than antennas seen in other species, possibly due to the low-light conditions in which 
these bacteria thrive. This antenna harvests sunlight, and then the excitation energy is transferred (via baseplate proteins, represented by the blue 
circles) to several FMO complexes. These act as a wire to transport the energy to reaction centers in the membrane. The reaction center shown 
here is from Chloroflexus aurantiacus ||73|| . The red hexagons represent the chlorosome envelope. Below this, in (b), there is a rendering of the 
structure of one of the FMO pigment protein complexes. The function of FMO in green sulfur bacteria is to transport excitations from the sunlight- 
harvesting LHC antenna (chlorosome) to the RC where it serves to initiate a charge-separation process. The FMO complex consists of eight 
bacteriochlorophyll a (BChl a) which are bound to a surrounding protein scaffolding. The excitation received from the antenna (here at molecule, 
or site, 1) is then passed from one BChl molecule to the next, until it reaches the molecule, or site, closest to the reaction center (here assumed to 
be site 3). It then irreversibly enters the reaction center, and ignites the charge separation process. The site labels here correspond to the states \j) 
is the Hamiltonian in Eq. Q and Eq. j2). 



T (1681]. Figure |2]3 shows the Liouville pathways, i.e. pulse-induced density-matrix dynamics of the system, that 
contribute to the beating signals. The time evolution of the coherence during the population time T has an oscillating 
phase factor, resulting in quantum beats in the 2D spectra as a function of T. The signal provides an incisive tool to 
probe quantumness of an electronic process. 

As we have described, the FMO complex 17711 is only one component of many types of light-harvesting complexes 
found in nature, but it is perhaps one of the most well studied and characterized. In the seminal work of Engel et al. 101, 
2D electronic spectra of the FMO BChl complex were obtained at 77K [see Fig.[3ja)], and strong quantum beating 
in the amplitudes and shapes of the diagonal peaks were revealed. These results have been interpreted ifisl [m 



as a strong signature of coherent wave-like transfer of excitation energy in the complex. More recently. 



experimental evidence has shown that quantum coherence survives in the FMO complex at physiological (ambient 
room) temperatures for at least 300 fs [7]. While much shorter than the coherence observed at cryogenic temperatures 
[see Fig. [3lb)], it is still thought to be sufficiently long to have an impact on the efficiency of the transport process 
(which, overall, is of the order of 1000 fs). This result is supported by a variety of theoretical models ll35ll which we 
will summarize shortly. Similarly, in the experiment performed by Collini et al. 101, two light-harvesting proteins, 
isolated from marine cryptophyte algae, were found to have relatively long-lasting excitation oscillations at room 
temperatures. This is a particularly interesting example, as the algae live in exceedingly low-light conditions, and 
must process the excitation energy transfer extremely efficiently to survive. In addition, the coupling of the LHC 
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Figure 2: (Color online) (a) Schematic of a 2DES setup, (b) Liouville pathways that contribute to the beating signals in the 2DES spectra. 
Contributions to 2DES signals are often depicted as double-sided Feynman diagrams as shown here. These diagrams represent the pulse-induced 
time evolution of the density matrix of the system. Light-matter interactions are depicted as arrows. An incoming arrow indicates an absorption by 
the system, whereas an outgoing arrow represents a stimulated emission event. Density-matrix elements are depicted in the center For example, 
in the R"'* pathway, the first pulse creates a coherence (superposition) between the ground state \g) and the excited-state \a), the second pulse then 
creates a coherence between \a) and ^S), and the third pulse prepares a coherence between \g) and ^8), which then emits the signal. The pathway 
will record a cross-peak at (w„, uip). Because during the population-time period (T) the system is in a excitonic superposition state, the signal will 
caiTy a phase factor that depends on T, leading to quantum beats on the 2D spectra. The rephasing and non-rephasing contributions are different in 
the time ordering of the first two pulses. 



systems to their protein environments differs in this case from most other photosynthetic LHCs, and Collini et al. 161] 
speculate that this may enhance correlated motions of the protein environment, enhancing any "coherence"-preserving 
effect that may arise. 

Because short pulses are used in 2D experiments, it is possible that the quantum beating signals in 2D spec- 
troscopy are results of vibration-coherence instead of electronic-coherence 11321 17811 . In principle, this alternative 
interpretation can be verified by examining the rephasing and non-rephasing spectra separately because vibrational- 
coherence would cause beating signals in all the rephasing and non-rephasing signals, whereas electronic-coherence 
would cause beating in the off-diagonal peaks of the rephasing signals or the diagonal peaks of the non-rephasing 
signals ll68ll . Recently, Turner and coworkers [|79l have carefully examined vibrational coherence effects and elec- 
tronic effects in 2D spectroscopy and confirmed that quantum beating observed in the PC645 light-harvesting antenna 
protein of the cryptophyte alga Chroomonas sp. is due to electronic coherence. 

In addition to 2D spectroscopy, other nonlinear optical experiments can be made into sensitive tools for quantum 
coherence phenomena. By using a two-color photon-echo approach, Lee et al. jstl also reveals that electronic coher- 
ence between bacteriopheophytin (H) and accessory bacteriochlorophyll (B) in the reaction center from the purple 



6 



bacterium R. sphaemides is preserved for much longer time than would be expected from the dephasing of either 
chromophores both at 77K and 180K Ql- Theoretical modeling revealed that the long-lasting coherence is indicative 
that the fluctuations of transition energies of B and H are strongly correlated in the reaction center It was speculated 
that the "coherence preservation effects" from correlated protein environments also explain the surprisingly long co- 
herence time (of the same order as the time for the excitation to traverse the whole complex) observed in FMO and 
other light-harvesting complexes. 



2.2. Theoretical models 

To gain a better understanding of what these various results imply it is helpful to consider a model commonly 
used to simulate some of these systems. One approach to simulate the quantum behavior of a single excitation in the 
FMO complex is via an effective Frenkel exciton model 13 [l^]- As described earlier, a photon creates an excitation 
in an antenna molecule, which is transferred to one of the sites in the FMO molecule, which is then transferred along 
the chain, from chromophore to chromophore via a transition -dipole Coulomb interaction, until reaching the reaction 
center, where it is then employed for charge separation. The assumption of a single-excitation seems valid for in 
vivo situations, particularly when the bacteria is living in very low-light conditions. It is exceedingly unlikely for the 
FMO complex to ever contain two excitations. Although the single-excitation assumption in experimental situations 
has generated certain controversy llssl . rapid exciton-exciton annihilation processes makes multiple-exciton states 
in a single complex extremely unstable and tend to relax into a single-exciton state before energy transfer occurs. 
Nevertheless, in principle it is straightforward to construct a multiple-excitation model, e.g., as employed in [80'] to 
describe an artificial photosynthetic system. For the purposes of clarity, however, we restrict ourselves here to the 
single-excitation model, which can be written as 



H^Yu + Z Jj.j'mj'\ + 

j=i i<j' 



(1) 



where the states I7) represent the presence of an excited electron (exciton) at site, or BChl molecule, j, where j e 1 , .., 7 
(see Fig. 1(b) for a figurative description of these sites in the FMO complex, and we omit the eighth BChl molecule 
because it was not discovered in structural models of the FMO complexes until recently). In other words, one can 
describe the whole FMO complex with a single label defining at which site the excitation is residing. The parameter 
€j is the energy of that excitation for a particular site (which is dependent on the surrounding protein structure and 
varies quite a lot), and Jjj is the excitonic coupling between the j and / sites. 

The determination of the various energies and coupling strengths in a given LHC is a compUcated and difficult area 
of research, often involving both spectroscopy and ab initio modeling of the physical structure and environment of 
the system In general, accurate values of electronic couplings that are in good agreements with experiments 

can be calculated with ab initio quantum chemistry methods based on the atomistic models of a LHC 1821 Isil 18311 . 
The site energies are more difficult to determine because modeling the protein-pigment interactions is a non-trivial 
task, therefore experimental inputs are usually required in order to obtain accurate energies. Recently, Renger and 
coworkers have demonstrated in several photosynthetic complexes 184 l85l 18611 that with a careful treatment of the 
electrostatic interactions between the electronic excitations and the surrounding protein environments, it is possible to 
calculate site-energy parameters with structure-based theoretical approaches. 

Since the FMO LHC is one of the most studied examples, some understanding of the system's Hamiltonian 
exists (though often different values for the energies and excitonic coupling amplitudes are given in the literature 
||87i|88n84J)- We now present the values used in Ref. 



H = 



215 


-104.1 


5.1 


-4.3 


4.7 


-15.1 


-7.8 ^ 


-104.1 


220 


32.6 


7.1 


5.4 


8.3 


0.8 


5.1 


32.6 





-46.8 


1.0 


-8.1 


5.1 


-4.3 


7.1 


-46.8 


125 


-70.7 


-14.7 


-61.5 


4.7 


5.4 


1.0 


-70.7 


450 


89.7 


-2.5 


-15.1 


8.3 


-8.1 


-14.7 


89.7 


330 


32.7 


-7.8 


0.8 


5.1 


-61.5 


-2.5 


32.7 


280 , 



(2) 



7 



where units are cm and the large energy gap relative to the common ground state has been subtracted. The diag- 
onal elements correspond to the energies e, in Eq. ([TJ, while the ofF-diagonals correspond to the site-site electronic 
couplings Jjj. We present these numbers here to show that the site-site couplings are of the same order as the energy 
difference between sites. Already this suggests that site-site coherences could be strong. In addition, the magnitudes 
of site-site couplings are often distributed broadly, with the larger ones in the scale of ~100 cm"', which is also 
comparable to the reorganization energy due to exciton-environment couplings. 



2.2.7. Approximations used in the models 

One of the open questions for this system is the nature of its interaction with its environment. In the simplest 
model, which we explicitly show below, a variety of approximations are made. The first is to assume that each site 
is coupled to a bath of oscillators, and that the baths coupled to each site are i ndep endent of each other The second 



approximation is to assume that this coupling is to the energy of each site ||69LI15LI901 |91L 19211 . which means that the 
environment causes fluctuations of the energy in the site basis. In the exciton basis, the fluctuations become to have off'- 
diagonal matrix elements that can cause transitions between eigenstates lfl6ll . The third is to assume the environment 
is Markovian (i.e., without memory). Finally, the fourth assumption is to assume the exciton-bath coupling is weak, 
and that one can perform a second-order perturbation theory to describe the dynamics of energy transfer 

This gives a master equation that describes the reduced density matrix of the FMO complex that is simple to solve, 
and has a minimal number of free parameters. Typically there is also a temperature-independent radiative relaxation 
rate for each site. If the excitation takes too long to traverse its way through the complex, it will be eventually lost 
due to this relaxation. Fortunately, fluorescence decay of chlorophylls typically occurs in the nanosecond time scale, 
a slow process compared to the other dynamics, ultimately giving the FMO complex its ability to efficiently transport 
energy with a high success rate. This simple model has estimated some of the qualitative properties of experiments 
on FMO, including the long coherence time, and high transport efficiency UtIi . 

As a specific example, the Markovian master equation approach assumes a self-energy, 

£[p] = Iss[p] + ^^dephM + S^dissM, (3) 

where p denotes the system density matrix. The first term describes irreversible excitation transport between site 3 



and the sink 1189119311 (see Fig. 1(a) and (b) for a description of where site 3 and the sink, or reaction center, are in the 
FMO): 

S^ssM = -^d[sds^i,P - 2s'^psD + PSds'd], (4) 

where sd - |3X0|, and To is the sink tunneling rate. Site 3 is coupled to the reaction center because it comprises the 
lowest energy excitation and lies closest to the reaction center. In contrast, site 1 lies closest to the antenna. Here jO) 
denotes the empty state with no electronic excitation. The second term, Sdeph[p], describes a temperature-dependent 



dephasing (see Ref. 119411 for a detailed discussion of the temperature dependence). The last term, Xdiss[p] describes 



the slow fluorescence relaxation process. A variety of insights have been obtained from this simple model 16^ 15 1 



which will discuss in the next section, but the essence is that the combination of coherent dynamics from H and 
level-broadening from 2deph[p] means there is a high possibility of bringing an excitation from site 1 to site 3, and 
thus the RC, before the excitation is lost due to the fluorescence relaxation Ediss[p]- The efficiency of the transport is 
often modelled via the time-dependent sink or reaction-center population 



Jo 



PRc(0=2ro ( Q''>(t')dt', (5) 

where Q^^^ - Tr[|3)(3|p(f)] is the population of the site "3", connected to the sink. 

As mentioned before, in most treatments of the FMO system, the coupling to the environment modulates the 
energy of each site, and in the language of quantum mechanics this is "pure-dephasing". The radiative relaxation 
of each site is independent of temperature because the optical transition energy of each site is exceptionally high 
(> 12, 000 cm '), leaving only the dephasing processes to be temperature dependent. 

In the FMO complex, each state in the above description represents an exciton on a molecule which transfers its 
energy to its neighbor, but no actual electron transport takes place. This implies that the re-organisation energy (the 
energy associated with a change in the surrounding protein environment due to the presence of an electron charge) is 
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Figure 3: (Color online) (a) The first example of quantum coherent beating between excitons 1 and 3 in the FMO complex at liquid-Nitrogen 
temperature, from Ref. j^l- (b) In another expeiiment (Ref. 01), signals of quantum beating demonstrate the temperature dependence of the 
coherence between excitons 1 and 3, from cryogenic up to ambient room temperature. Note that by inspection it appears that coherent beating is 
observable in the room temperature data for approximately 300 fs. (c) An example of the 2D electronic spectra of the FMO used to collate figure 
(b). The data for figure (b) is extracted from the point indicated by the white circle, which is expected to contain the inforaiation about coherence 
between excitons 1 and 3. The two-dimensional electronic spectra data is collected by applying a seiies of three pulses to the sample, and observing 
the output signal. Theoretical models indicate that the cross-peak (which an off-diagonal peak in terms of the whole spectral picture) in the signal 
spectra directly corresponds to the presence of coherence between different excitonic sites. Mathematically this is equivalent to a kind of quantum 
process tomography fTOlpTlll , or a probing of the memory of the system. 
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Figure 4: (Color online) This figure, from Wilde et al BtIi . shows a theoretical violation of the Leggett-Garg inequality as a function of dephasing 
(and hence bath temperature) for measurements on each site, 1 to 7. A value of < indicates a violation, and the results suggest quantum 
phenomena persists up-to room temperature for measurements on sites 1, 2, 5, and 6. As yet there is no way to directly measure such non-invasive 
observables in FMO. It is particularly difficult, as unambiguous violations of the Leggett-Garg inequality in this case would require non-invasive 
measurements of excitonic populations. Most fluorescence and spectroscopic measurements are destructive, and thus not appropriate for observing 
the Leggett-Garg inequaUty. However, ultimately this kind of test would help rule out alternative equivalent explanations for efficient transport in 
FMO. 



relatively low compared to processes where electron transport occurs. In general, however, both the assumption of 
independent Markovian baths and of weak coupling to these baths may not be realistic. There is evidence that the 
coupling to protein environments around the FMO complex can be of the same order as the electronic couplings ll94l 
[sol] (100 cm"^'). In addition, the bath (e.g., the protein scaffold surrounding the FMO complex) may have structure 
and dynamics which have a strong correlation with, and back-action on, the dynamics of the excitation in the FMO 
complex. To understand the effect of this intermediate coupling regime, and complex environment, the system and 



bath have been treated with a range of non-Markovian and higher-order (in electron-phonon coupling) models 1135 
0, Issi In particular, nonperturbative approaches that provide exact numerical results for 

certain models o f exc i tation energy tr ansfer have been applied to study coherent quantum dynamics in photosynthetic 
light harvesting I lOOi 101, 102 , 103 1. We will discuss the details of one of these approaches (the hierarchy model) 
later Although these nonperturbative methods are computationally too expensive to apply to large photosynthetic 
complexes, they have provided valuable insight that have shed a different light on the coherence effects in light 
harvesting in a different light ll35ll : one where bath memory effects conspire to enhance quantum coherence. 

Going beyond the independent bath approximation, initial experiments indicated that the surprising long-lasting 
quantum coherence between two electronic states may be because of coherence enhancing effects from coupling 
to commo n vibrational modes |69, 5lB| - Thi s has triggered the reexamination of the excitation transfer processes 



m 



10711 in the presence of such an unusual environment. In princi ple, s uch a 

Iioslliogli . 



1 104 110511 and decoherence effects 

strongly coupled "common bath" should contain bath-induced fluctuations that are spatially correlated 



As a c ounte r-argument, however, some molecular dynamics simulations for the FMO complex II34I1 and reaction 
center show that only weak correlations between the movements (vibrations) of the chromophores appear, and 
that the uncorrected bath approxim ation s may be valid. 

Very recent work by Shim et al ll 1 ill compared some of these master-equation models to a more complex many- 
body atomistic/molecular modeling scheme and found excellent agreement between the two approaches. This is 
encouraging as it does suggest that the assumption (based on relatively simple models) that the oscillations seen 
in experiment have a truly quantum origin might be correct. However, whether the environment surrounding LHC 
systems has a structure or nature that preserves quantum effects in some way is still a controversial and open question. 
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More experimental work is required to answer this puzzle. 

Finally, some authors have argued that it is p ossible to construct a variety of alternative classical models that can, 



in principle, produce both classical beating 113311 and efficient energy transport 11691 11511 without reference to quantum 
coherence at all. Such alternative descriptions must be eliminated before one can unambiguously state that the FMO 
complex, or other LHCs, take advantage of quantum mechanics. 

2.3. New Insights from Theoretical Modelings 
2.3.1. Environment-assisted transport 

Apart from the direct observation in experiments of quantum coherent beating, some of the most intriguing in- 
sights that have appeared from the study of the FMO complex suggest that a combination of these quantum coherent 
oscillations of excitations between sites, and interactions with the environment, produce a transport efficiency higher 
than is possible with a Forster model [0] alone. If correct, this would fit the definition of functional quantum biology. 
Such a phenomenon was proposed and studied by Plenio et al lllSll . by Mohseni et al andbyLeeefaZ iH. In 



all cases, it is assumed that a single excitation is placed at a particular site in the FMO (usually site one), and that this 
excitation then propagates through the FMO chain due to a combination of coherent tunnelling and environmental 
effects. As mentioned earlier, the goal of the FMO complex is to get this excitation efficiently to site three, which 
is coupled (incoherently) to a reaction center Thus the typical time, or efficiency, of reaching the reaction center is 
calculated as a function of environment temperature, site selective-couplings, and so on. In most cases it seems that 
the full quantum model (e.g., based on Eqs. [1-3]) will give a higher efficiency than a purely classical one. Several 
physically intuitive explanations have been given for this phenomenon: 

• Local minima avoidance: Refs. | f5] and lfl6ll showed that the Markovian quantum model outlined earlier sug- 



gests both long-lived coherent oscillations between sites, and an enhanced rate of excitation transfer to reach 



the "sink" site, over that predicted by the Forster model. Reference 111511 argued that essentially the combination 
of coherent transfer and environmental "noise" causes level broadening, which implies that the excitation can 
more easily escape local minima in the FMO network (see Fig. [ 1 ]). In other words, quantum delocalization can 
help avoid and overcome local minima, or energy potential traps, in the energy landscape of the FMO complex 
13511 . It has been suggested that such a phenomena is even more important in higher plants beca use o f possible 



"up-hill" potential landscapes. These results are related to earlier work by Gaab and Bardeen Hi 1211 who out- 
lined how energy transfer through a network can be optimized by carefully choosing the coherence time and the 
rate of transfer to the sink. 

• Coherence assisted trapping: In a study by Lee et al. lfl9ll . it was demonstrated that coherent quantum dynamics, 
together with rapid incoherent dissipation due to a trapping site can enhance the efficiency of the irreversible 
energy transfer between the initial energy donor and the sink, or reaction center, over a purely classical model. 
Effectively, this model suggests that an intricate combination of reversible quantum coherent evolution and 
incoherent collapsing of exciton wave functions can achieve a greatly enhanced energy-trapping efficiency, 
effectively using the anti-Zeno effect to promote energy transfer 

However, one should again notice that this "enhanced efficiency" is in comparison to that obtained with the Forster 
model. There remain several ambiguities; Firstly, is the Forster model the correct classical model to make a compar- 
ison to? Secondly, placing a single excitation at site "1" may match well certain experiments, but in normal light 
conditions is this a correct approach? 

Quite recently, Briggs and Eisfeld ll33ll showed that, for realistic coupling strengths, an alternative classical model 
of the FMO transport process can produce results identical to the quantum one. In addition, a recent analysis by Wu 



et al. Ill 1311 suggests that, as mentioned earlier, the efficiency enhancement gained by coherence in the FMO complex 
is only a few percent compared even to that predicted by the Forster model. Finally, the true in-vivo conditions are 
not well understood, and is the subject of current research. Perhaps the only way to really sol ve this issu e is to show 



that in more complex components of some light harvesting complex, e.g., in LHl or LH2 111 14 Il09tl . there exist 
significant energy traps which, without the assistance of quantum coherence, drastically impact/reduce the probability 
of an excitation successfully navigating its way to a reaction center before being lost to fluorescence relaxation. 
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2.3.2. The Hierarchy model 

In an effort to gain a deeper understanding of the interplay between the quantum coherent transport of energy in 
FMO and how it interacts with its complex protein environment, a variety of non-Markovian and non-perturbative 
models have been appli ed to these systems. One of the most successful is the Hierarchy model, originally devel oped 



by T animura and Kubo 111 1511 . which has been applied to both the FMO complex 0511 and other LHC components 111 14 



109ll . This model has a had a large impact this field, and thus w e will briefly describe its main components here. Full 



derivations can be found in the literature 111 16[ ll 17l ll 15l ll 18ll . Starting with the Hamiltonian we discussed earlier, 



one explicitly describes the interaction term between site energies and phonon modes, 

H^''-"' = 2 1;)01 qj (6) 
j 

where qj - -Y^kgj.kXk, and gj^k is the coupling constant of the jxh site and the k\h mode in the bath. The phonon 
modes themselves are described as harmonic oscillators (with a Hamiltonian H*^''' which we omit here for brevity). 

It is typically assumed that at f = the sites (or pigments) and the phonon modes are separable, so that p(0) - 
Pe(0) ® P/,(0), and the phonon modes are in a thermal equilibrium state pp(0) = e-^«""/Tr[e-^«""],j6 = IIK^T. After 
averaging, the correlation function of the bath modes 

1 r°° e~"^' 

Cj{t)^{qj{t)qjm^ - dojJ,Au)- 5- (7) 

71 Jo I — e p 

are sufficient to describe the properties of the bath. With the Hierarchy method one typically uses a Drude spectral 
density (appropriate for over-damped oscillators) 

Here is the "Drude decay constant", and indicates the memory time of the bath for site j (each site is assumed 
to have its own i ndep endent bath, though in general one of the powers of the Hierarchy method is its ability to treat 



correlated baths 110911 ). Also, Aj is the reorganisation energy, related to the system-bath coupling strength. This Drude 



spectral density implies an exponentially-decaying correlation function, 

= Yj ^i'"' ("^ i'"'') 

«i=0 

where - jj, fij,„ > 1 - Ijimlhfi, and the coefficients 

C;,0 = Jj'^j {cot(fihyj/2) - i) /h (10) 



and 



With just this information on the bath properties (the reorganisaiton energy, bath memory time, and bath tem- 
perature) one can employ the Hierarchy equations to describe the system-bath dynamics in the strong-coupling and 
non-Markovian regime, 

N K N K 

p„ = -(/L + ^ ^ nj,mM^m)Pn ''YjYj t^-''^"!"] 

j=\ ff?=0 j=l m=0 

N K 

- ' X X ('^'" ~ ^j) ■ ^ 

j=l m=0 
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Here, Qj - \j}{j\ is the projector on the site j, L is the Liouvillian described by the Hamiltonian introduced earlier, 
and for FMO N - 1 . The various parameters are those we defined above. The Hierarchy is a large set of coupled 
equations each labelled by n, a set of non-negative integers uniquely specifying each equation. The integers are 
defined as n = {111,112, n^, 11^] = {{«io, «ii, •■, nix], ■■, [nm, nN\, iink]]- That is, each site j has an additional label 
m, from to K, and each of those labels in turn can run from to 00. The label n = = {{0, 0, 0....)j is special, and 
refers to the system density matrix. Its properties at any time t define those of the system. This is in turn coupled to so 
called "auxiliary density matrices", which describe the complex bath fluctuations, by the terms in the equation where 
n^ijj (i.e. n^jjj implies the term in the index defined by j and m is increased or decreased by 1). As a simple example, 
if the terms in the sum for m are truncated at K - 0, then the system density matrix is coupled to other auxiliary 
density matrices directly, Oj = {1,0,0, ...0), O2 = {0, 1,0, ...0), and so on (see Fig. |5jc)). Each of these auxiliary 
equations is then in turn coupled to higher-tier equations, and overall the hierarchy forms a high dimensional simplex. 

The hierarchy equations are infinite both in the sum of the m label (due to the expansion of the correlation func- 
tion), and in the value of each label itself. Thus both must be truncated in some way. Typically, the m label can be 
truncated if jShy,,, < 1 and 2n(K + l)/fifi > cos (where ais is some characteristic system frequency). The overall 
truncation of the labels is defined by the largest total number of terms in a label, or the tier, Nc - Yji,m nj,m- This 
conc ept is important as each ADO in a given tier is only coupled to those ADOS in the tier above or below it. In Yan 



et al 111 1811 they describe this structure, and the meaning of each tier, as the level of bath self-correlation included in 
the simulation. 

The tier at which the truncation should be made is difficult to predict (there is no small parameter here), though 
typically truncation should be taken when Nc-y » cos (though it also depends on the reorganisation energy and tem- 
perature). In practice, truncation should be checked by looking for convergence in system dyna mics when changing 
Nc- A systematic and powerful algorithm to find convergence was introduced by Shi et al jl 1911 . They proposed that 
a renormalization of the hierarchy equations allows a direct inspection of the elements of higher auxiliary equations 
to be made during the numerical solution, and thus the level of the tier can be truncated as necessary. 

In addition, a temperature correction can be added to the equations of motion, which compensates for the trun- 
cation of the m label at level K, and assists in overall convergence. This is made by assuming that all higher-level 
auxiliar y den sity operators for m > K are uncoupled from each other and that they experience pure Markovian dy- 
namics il 16ll . This r esults in an additional term in the equation of motion, sometimes called the "Ishizaki-Tanimura 
boundary condition" Il 16ll . which can be included in the Hierarchy equation of motion; 

A' 00 

^-Bc = -2 Z ^[e.ie.Pn]]. (13) 

j=l m=K+l ^J'"' 

The double commutator gives a normal Lindblad form, and the summation can be done analytically, e.g., for K = Q, 

y — = -^^[1 - yjh[coi{yjhlill2))lil2]. (14) 

In other words, this term encapsulates additional environmental dephasing effects neglected by truncating K. Including 
this correction, even for K - Q, has been found to substantially assist the convergence of the Hierarchy results. As 
an example, in Fig. (|5]i, we show the results of an application of this model to the FMO system at 77 K and 300 K, 
for various yj and Aj employed in the literature. The beauty of this method is that, to some degree, all parameters can 
be extracted from experiments and the results correctly predict the long coherence time seen in experiments. Perhaps 
the most important point here is that long bath memory (implied by small yj) helps increase the coherence time of 
the coherent oscillations llssll . though this may not necessarily enhance the overall transport efficiency. In general, for 
each parameter, there is an opti mal value which gives the largest transport efficiency from site 1 to site 3, and hence 
to the reaction center 



2.3.3. Leggett-Garg inequality and entanglement measures 

As mentioned earlier, it is important to fully verify whether the experimental observations of coherence truly 



come from quantum mechanics and not some alternative, classical, model [33]. Wilde et al 13711 proposed using the 



Leggett-Garg inequality as a means to unequivocally ascertain whether the dynamics observed in these experiments 
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Figure 5: (Color online) (a) The population dynamics predicted for the FMO complex by the hierarchy model. The parameters used were =50 
fs and /i = 35 cm" 1, for all sites, as in (S^l . The dashed lines indicate a bath temperature of 300 Kelvin and solid lines 77 Kelvin. Truncation was 
taken at K = and N,. = 4. The lower figure is a representation of the sparse matrix of coefficients describing the coupled hierarchy equations. For 
Wj. = 4 and N = 7, there are (N^ + Ny./(Nr\N\) = 330 total density operators in the hierarchy. Each density operator obeys the dynamics of the 
7-site model, and thus in a typical open-quantum-system super-operator basis is described by 49 equations. In total this gives a matrix of equations 
of size 16, 170. It is helpful to present this matrix in (b) and (c) as it serves as a visual check for faults in implementation (e.g., due to missing 
symmetry or so forth). The magnification of the 8 X 49 = 392 elements, in (c), illustrates how the system density matrix is coupled to the seven 
A'c = 1 tier of auxiliary density matrices (the coupling terms are highlighted in black). 



are quantum or obey "macroscopic realism". As yet, however, the measurements needed for such a test cannot be 
performed in these FMO experiments. However, their results show that if it were possible to perform such measure- 
ments, the currently accepted physical parameters for the FMO system suggest a violation of the inequality at room 
temperatures should occur (see Fig. 3} for an example of the violation of the inequality. 

In addition to coherence, it is also interesting to ask whether there is ap prec iable entanglement in the FMO complex 



during the exciton transfer It was recently suggested by Sarovar et al. Il36[] that a small amount of long-range and 



multipartite entanglement should exist even at physi olog ical temperatures 113611 . however, the role of entanglement in 



the energy transfer process is still not clear (see also [|31|). 



2.3.4. Quantum process tomography 

Quantum state tomography is a well-known technique to reconstruct the full density matrix or state o f a sys- 
tem. Recently, there has been increasing attention on what is termed quantum process tomography llM [ml, which 



attempts to fully characterize a given quantum process (i.e., a systems dynamics) in an open environment (e.g., in 
contact with Markovian and non-Markovian baths). Simply speaking, quantum process tomography is a systematic 
procedure to characterize an unknown quantum system, or black box, by systematically analyzing the functional rela- 
tionships between inputs and outputs as a function of time. Very recent works [T^ 71 1 showed that the 2D electronic 



spectra used to investigate light-harvesting complexes is equivalent to a kind of quantum process tomography, eluci- 
dating the meaning of the results found with these techniques in earlier experiments. In principle this means that the 
spectroscopic experiments of excitonic systems can be studied from a quantum information theory approach, and that 
perhaps the full nature of the dynamics in FMO and other systems can be understood on a deeper level. In the future 
perhaps it may be possible to discuss important questions such as the existence of non-Markovian environments, en- 
tanglement, and decoherence during the transfer process using this tool. However, their proposal still requires steps 
that are experimentally challenging. 

2.4. Functionality 

One of the arguments used to justify the presence of "functional quantum mechanics" is that if there is any 
advantage, or improvement in efficiency, to be had from taking advantage of quantum mechanics, evolution will 
inevitably do so. Quantum coherence, manifested in the delocalized eigenstates of photoexcitations in photosynthetic 
complexes, plays fundamental roles in the spectral properties and energy-transfer dynamics of photosynthetic light 
harvesting |[l2i [Till . In structures of strongly-coupled clusters of pigments, coherence is essential for shaping the 



14 



energy landscapes and energy flow tow ards the reaction center Thi s has been demonstrated clearl y in the LH2 
antenna complex of purple bacteria |59, 122 , 65 1, t he F MO complex 1 123 , 89ll , the reaction centers 1 124 1. and the 
major light-harvesting complex of higher plants fesl Il25ll . For example, in the LH2 antenna complex, energy transfer 
from a B800 BC hl to coherently delocalized B850 states is an order of magnitude faster than that to a single B850 



BChl 11124 112211 . Such a significant enhancement in energy transfer rate clearly contributes to the highly efficient 



light-harvesting process, hence the functionality, in purple bacteria. 

A separate issue is whether or not quantum coherent dynamics, i.e. excitonic coherence, plays a role in the func- 
tionality of photosynthesis. This issue has become the subject of intensive research rec ently. How ever, in reality the 
calculated gain in efficiency (as a function of time, for example) seen in Refs. |fl5 [ [Tel. [l9i [l0i [TTsIl is relatively small. 
In some species (e.g., the algae living in low-light conditions studied by CoUini et al j^) this small advantage may 
be useful. But in general, (e.g., plants in high-level light conditions), is there still a significant advantage provided by 
coherence? As mentioned earUer, some studies of "higher plants" suggest the energy landscape of their LHCs is even 
more rugged than that in the FMO complex fl, and thus coherence may become even more important in efficiently 
navigating this terrain before other processes (like fluorescence relaxation) steal away the excitation. Interestingly, a 
recent study by Hoyer et al. U2m indicates that the effect of "quantum speedup" is rather small in the time scales of 
photosynthetic energy transfer, and they propose that quantum coherent effects in light-harvesting systems are more 
likely to be optimized for robustness or total efficiency. 

Finally, very recent results by Ritschel et al lISSll paint an even more complex picture; they model the full FMO 
trimer, including an additional eighth bacteriochlorophyll (BChl) molecule which was recently found to exist in the 
FMO complex. On a positive note, they find that the subunits of the trimer act as separate transport channels from 
the antenna to the reaction center, thus the assumption of studying the individual FMO monomers individually seems 
safe. However, if they use this new eighth site as the source of initial excitation (which may be the case, as it is located 
close to the antenna), they find that no coherent exciton dynamics occur; only exponential decay. In addition, they find 
discrepancies in the overall transport efficiency, depending on which energies and couplings (i.e., which Hamiltonian) 
are chosen from the literature. In contrast, Olbrich et al. Kill used a Molecular Dynamics approach to derive a time- 
dependent equation of motion for the excitation transfer, and found that the eighth BChl did not form a stable complex 
with the other seven, and that, while thermal effects in total were stronger than assumed in other studies, the coherent 
effects were enhanced when considering the trimer system in totality. This is attributed to a corresponding increase of 
the site coupling within a monomer due to the influence of the other monomers in the trimer However, they found that 
the quantum coherent effects were still quite small in their model, and Olbrich et al. 1 8 Til argue that the fast thermal 
fluctuations of the site energies are the main contributing factor to the efficient transport of energy. 

Inevitably a broader picture is needed before we can fully understand the role of quantum coherence in natural 
conditions. Is the coherence observed in experiments contributing significantly to the efficient transfer of energy in 
natural conditions? In-vitro experiments and theoretical models suggest a positive answer, but ultimately experiments 
into other types and classes of LHC systems (e.g., LH-1 and LH-2) are required, and a broader understanding of the 
transport process in vivo is needed. 



3. Avian Magnetoreception: a tale of two spins 

Magnetoreception is a unique animal ability to detect either the polarity or inclination of the earths magnetic 
field as a navigation tool. As with electro-reception (the ability to detect electrical fields, e.g., sharks hunting their 
prey), it represents a sense slightly alien to human experience. Or so it was thought, until recently. Experiments 
on manipulating a specific gene in flies and monarch butterflies have indicated this gene plays a significant role in 
magnetoreception, and since humans share this gene in a dormant form it suggests that even we may harbor this 
unique ability ll52i Issll . However, some evidence [Issll suggests it is not driven by the radical -pair mechanism we 
discuss here. 

The epitomic example of magnetoreception however is in migrating birds. They apparently use this powerful abil- 
ity, alongside other biological geo-mapping tools, to navigate large distances. A complete and satisfying description 
of the mechanism for the magnetoreception ability of many species remains a mystery. Some evidence suggests that 



a magnetite-based (deposits of magnetic iron minerals) mechanism describes well this ability in some animals 143 1 



(though recent r eport s suggest that, at least in pigeons, these deposits are in fact macrophages and play no role in the 



magnetic sense Ill27ll ). In other avian species, particularly robins, several intriguing aspects of their behavior rule out 
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magnetite-based mechanisms. Here we focus on the radical-pair model for magnetoreception that represents another 
potential functional quantum biological system. Precise details will be given later, but summary radical pair is a pair 
of molecules bound together, with each has an unpaired electron. These electrons form naturally into singlet or triplet 
states, the state of which, in turn, affects the disassociation, or reaction products, of these radical pairs. It is proposed 
that the relative weights of these reaction products which are biologically detectable. Thus, if the relative probabilities 
of having singlet or triplet pairs is effected by the Earth's magnetic field, this in turn affects the distributions of the two 
different reaction products. As a result, a biologically detectable signal can be generated to function as a compass. 



3.1. Behavioral studies 

The hypothesis thatmigrating birds can make use of the earth's magnetic field for orientation was first pro pose d 
by von Middendorff lll28ll . Later, European robins were found to utilize a magnetic compass when migrating lll29|] . 
Since then, many experimental tests on the behavior of migrating birds have been performed. In particular, it was 
found that for some species their magnetic sense acts as an inclination compass, but is insensitive to parity H. In 
other words, they can detect the relative angle of the local earth's magnetic field lines, but not whether they are going 
north or south. There is some benefit to this in that the compass would not be effected by geomagnetic reversal lll29ll . 
Typically, however, most models of magnetite-based magnetoreception act as a parity compass ||43|. 

With advances in technology and technique, many more delicate behaviorial tests have been performed. For 
example, it was found that the compass was: (a) dependent on the frequency and intensity of ambient light fl42l 14311 : 
(b) dependent on the intensity of the external magnetic field ll43ll : (c) disrupted when repeated magnetic pulses were 
applied to the migratory birds ll44ll . and (d) disrupted when an external oscillating field was applied at certain angles 
to the geomagnetic field 14511 . In additi on, la teralization in terms of the the birds' visual system associated with 
magnetoreception has also been reported 111301, suggesting the mechanism operates primarily in the right eye. 

All of this data lends support to a photo-activated radical-pair mechanism. In particular, the experiments showing 
(d), the disruption of the navigation mechanism by an external oscillating field in the megahertz range (as seen in 



migratory robins and directionally-trained chickens as well as zebra finches 11431 145L 1471 14611 1 match exactly what is 
predicted by simple models of radical-pairs, which we will describe shortly. 

Very recent evidence (see Ref. ll43ll for a detailed review) suggests that both an inclination compass and a polarity 
compass are present in certain avian species, suggesting that a magnetite and radical-pair compass can coexist. We 
also point out, as mentioned earlier, that in some species of flies, butterflies, and even perhaps humans, there may be a 
third unknown mechanism at work [S^ 53] that has some of the same behavioral properties of the radical-pair model, 
but not others. Here we will focus on describing the evidence for the radical-pair mechanism, and in what way that 
mechanism represents "functional quantum biology". 



3.2. Radical-pair model: evidence and motivation 

Based on findings of behavioral studies, it is obvious that a successful model should contain a light-dependent 
factor Another feature is that because the geomagnetic field is very weak (50 /^T), the magnetic compass mechanism 
must be very sensitive. In addition, this mechanism must transduce the response to the Earth's field into a biologically- 
detectable signal. Finally, as we mentioned earlier, in some species the compass is only an inclination compass. 

To help satisfy some of these criteria, the radical-pair mechanism (RPM) for magnetoreception was proposed by 
Sch ulten et al. 13911 in 1978. Within the radical-pair model there exists the possibility of a photo-sensitive element 
1 4llll3lll . i.e., absorption of light triggers an electron transfer from a donor to an acceptor molecule to form the initial 
radical pair (biologically this is proposed to occur within the cryptochrome proteins within the eye). This process can 
be made to satisfy most of the frequency and intensity dependencies observed in experiments on Robins 1143 1. 

How then is this donor-accepter system sensitive to the Earth's field? Most known radicals contain atoms (hy- 
drogen and nitrogen) with nuclear spins which act like an effective internal magnetic field affecting the electron spin 
via the nuclear hyperfine interaction. To sense the weak magnetic field of the Earth, the combination of anisotropic 
internal nuclear hyperfine fields and external geomagnetic fields causes a mixing of the electron singlet and triplet 
states, the strength of which is dependent on the angle of the external geomagnetic field. Ritz et al lll3lll assumed 
a simplistic internal magnetic environment with only one anisotropic nucleus (i.e., one radical is devoid of internal 
magnetic fields, and the other radical has very stro ng in t erna l magnetic fields) and showed that this maximizes the 



singlet-triplet mixing compared with other designs 11311 113211 . Whether this occurs in nature is unknown, but their 



extreme example helps to present a simple model of how this mechanism works (see later). 
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Figure 6: (Color online) (a) The radical-pair reaction mechanism as proposed for the avian compass. There are three main steps in the usual process 
of magnetoreception via this mechanism: (1) Light-induced electron transfer from one radical pair-forming molecule (e.g., in the cryptochrome 
in the retina of a bird) to an acceptor molecule creates a radical pair (2) Singlet (S) and triplet (T) electron-spin states inter-convert due to the 
the external (Zeeman) and internal (hyperfine) magnetic couplings. (3) Singlet and triplet radical pairs recombine into singlet and triplet products 
respectively, which are biologically detectable, (b) Singlet yield as a function of external field angle 9 in the presence of an oscillatory field, but 
no noise (from |23|| ). The blue curve shows just a static geomagnetic field (Bq = 47 /jT), and the red curves show the singlet yield in the case 
where a 150 nT field oscillating at 1.316 MHz is superimposed perpendicular to the direction of the static field. An appreciable eft'ect on the singlet 
yield occurs once k (the decay rate of the radical, equivalent to as defined in equation 1 17) is of order 10"'' .s"' . (c) Singlet yield as a function 
of magnetic field angle 6 for differing noise mag nitudes (from HI) (for k = 10*). The blue curve shows the optimal case with no-noise. The red 
curves indicate that a noise rate of T > O.Ik has a detrimental effect on the singlet yield. Both of these results indicate that the electron spin state 
must have a remarkably long coherence (or congelation) time. 
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For t he last stage, as mentioned, this radical pair can decay into different reaction products at different rates 



I 41LI133I1 depending on the spin state (singlet or triplet). If this radical-pair is located in the cryptochrome proteins 
within the eye, this could lead to a vision-based compass where the bird actually sees a directional signal [14 lH . Finally, 
this model requires that the many molecules forming the rad ical-pair-b ased magnetic compass must be ordered in some 



u-p: 

way so that the directional effects will not be averaged out II 134 1 13511 . This may require that the radical-pairs exist in 



a regular ordered biological structure, or lattice, which has not yet been identified. 
3.2.1. The Schulten model 

A simple model to calculate the influence of the Earth's field on the reaction products of the radical pair was 



outlined originally by Schulten 09114111 . Here we present specific details of this model, to help clarify the physical 
picture described in the previous section. As discussed above, the radical pair is formed when an electron is transferred 
from a donor to an acceptor molecule, typically via an optically-excited process. This forms a singlet state for the 
combined donor-acceptor system. The, ideally the subsequent dynamics of this state is primarily affected by an 
anisotropic hyperfine interaction between the electron spin and the nuclear spin within the molecules. This must 
occur faster than any decay (recombination) or decoherence of the radical-pair The interaction is described by a 
hyperfine coupling Hamiltonian (assu ming only one of the two spins is affected by its internal nuclear spin, for the 



most sensitive possible compass 1)13111 ) 

Hhyp^Si-A-I, (15) 

where S i - (cr^/2, o-y/2, (t^/2) is the operator of electron spin 1, / is the nuclear spin operator, and A - (A^, A,,, A,) 
in its diagonal basis is the anisotropic hyperfine coupling tensor for a cigar-shaped molecule with A , - Ay - A, 12 
and A, = 2.418 MHz, for high sensitivity lEsIl (though recent results have suggeste d th e weak-hyperfine coupling 
regime can produce a stronger sensitivity, given the right combination of parameters lll36il ). The Zeeman splitting of 
the electron spins is 

Hn^gtiB^-iSi+S.), (16) 

where B is the magnetic field vector describing the external geomagnetic field (or any other applied fields), g - 2, 
and is the Bohr magneton of the electron. The magnitude of the geomagnetic field is typically of the order 
Bo ~ 0.5 G (~ 50 juT). For the axial symmetry of the hyperfine tensor, the magnetic field vector is explicitly written 
as B = Bo(sin6',O,cos0), for 9 e [0,;t/2]. 

In addition to the coherent evolution driven by Hhyp and Hb, the radical-pair dynamics feels several incoherent 
effects. Both singlet {S) and triplet (r+. To, r_) states undergo spin-selective relaxation or recombination into the 
singlet product state |s) and the triplet product state |t). This occurs when, e.g., there is electron back-transfer from 
the acceptor radical to the donor radical, which is described by the self-energy 

^^plp] = -Tr ^[SapslpP - 2slppSal3 + pSapS^p], (17) 

where sip - \S,/3} (s|, ^2/3 = \T+,/3} (t|, ^3/? = \To,/3} (t|, and 54/3 = \T-,/3} {t\, for the up and down nuclear spin states 
/3 =t- i- Here we assume that the singlet and triplet products have the same recombination rates Fr, = Fr, - Fr. The 
final incoherent term is the environmental noise. The simplest noise model is where both phase and amplitude of both 
electron spins are equally perturbed with self-energy 

^n[p] = -r^YjlLkjLljP - 2Ll.pLkj + pLkjLl], (18) 

where Lij = cr ^, Lij - cr,,, - cr- for the electron spins j - 1,2. 

The dynamics of the electron and nuclear spins are then determined by the Hamiltonian H - Hhxp + Hb together 
with the incoherent contributions S = 2s+2nforthe generalized master equation: p(f) = X[p(0]> where the Liouvillian 

i; is X = -i[H, ■] + s. 

The biologically-detectable signal is the integrated triplet or singlet yield. For example, the triplet yield is, 



Jo 



I rR,T{t)dt, (19) 
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where, T(t) = Ti[Q^pit)], and is the projector onto the triplet states. The interplay between the external geomag- 
netic field and the anisotropic nuclear hyperfine interaction results in a dependence of the triplet yield on the relative 
angle between the geomagnetic field and the radical-pair molecule. 

This simple model is suflicient to describe the destructive effect of additional radio frequency fields, as observed 
in experiments, and enables one to estimate the most ideal properties of the radical-pair system to maximize the 
sensitivity of the magnetic compass. 



3.2.2. Radical-pair molecules: experimental tests 

As discussed earlier, a variety of behaviorial experiments support the radical-pair mechanism; e.g., the disorien- 
tation of the magnetoreception ability in Robins when external magnetic fields of around a megahertz frequency are 
applied 137 , 471. In addition, comparison of the frequency-dependence of this behavior to theoretical models 
indicates that the radical-pair system does have the optimal design mentioned above; i.e., one of the radicals may 
contain almost no nuclear spin polarization. 

To further bolster evidence supporting the radical-pair model, there are also experimental tests at the molecular 
level on a variety of candidate molecules that could be at the heart of the mechanism. There have been a large amount 
of studies on radical-pairs in solution, and on the efifect of external fields on the reaction products (see a review of 
early experiments by Steiner and Ulrich IHSll ) though only under magnetic field strengths much stronger than Earth's. 
In recent work. Woodward et al ll49ll developed a technique for studying the properties of such candidate reaction-pair 
systems. As an example they studied the photochemical reaction of pyrene (Py) and N,N-dimethylaniline (DMA) in 
solution under a very weak (500 /iT) radio-frequency magnetic field. When pyrene is continuously irradiated by UV 
light, in the presence of DMA, the transient radical ion pair Py^ DMA^ is produced in a spin-correlated singlet state 
and is found to interconvert with the corresponding triplet state, as expected. An exciplex (an electronically-excited 
complex of Py'" and DMA' ) is then produced due to the electron-hole recombination of the singlet radical pair Their 
technique uses the fluorescence of the exciplex to monitor the concentration of singlet pairs, and estimate the reaction 
products. They found that even this very weak magnetic field induced changes of up 25% in the reaction products. 

Rodgers et al. ifsoll also studied the spin-selective recombination of Py " and DMA " radicals under the affect of 
0-23 mT magnetic fields. It was found that weak magnetic field effects become most pronounced when the ratio of 
hyperfine coupling strengths for the two radicals is large. Another recent example is the triad composed of linked 
carotenoid (C), porphyrin (P) and fullerene (F) groups considered by Maeda et al. 15111 . By using spectroscopic 
measurements they found this particular radical-pair's lifetime was changed by the application of magnetic fields less 
than 50 yuT (as needed to detect the Earth's weak geomagnetic field). In addition, they found this eff'ect responded 
anisotropic ally to the angle 9 between the molecule and the external field, which is of course an essential feature of a 
radical-pair compass. However the anisotropy was only manifest at stronger fields (3.1 mT), and no anisotropy could 
be observed for 50 /iT, even at low temperatures (113 K). It is hoped that future work will show whether a radical 
pair formed by "reduced flavin cofactor and an oxidized tryptophan residue in a cryptochrome flavoprotein" (Maeda 
et al ifsill ) will have the long spin coherence lifetime and low recombination time to achieve the 50 /yT sensitivity at 
room temperature needed for this mechanism to function. 



3.3. Functionality and significance of quantum coherence 

Attention has recently focused on the functionality of quantum coherence in the radical-pair model because, as 
we noted in the introduction, spin singlet and triplet states are inherently quantum. This has led to a great degree of 
interest in this phenomena from researchers in condensed matter and quantum information theory. For example, an 
analysis of the entanglement and decoherence effects of this mechanism 112 311 suggests that, in their words, for this 
model to function as desired, ''superposition and entanglement are sustained in this living system for at least tens of 
microseconds, exceeding the durations achieved in the best comparable man-made molecular systems" 

This implies that the radical-pair model is, if it really exists and functions as described, a piece of functional 
quantum biological hardware. From the perspective of condensed matter physicists, it is very akin to man-made 
artificial molecule systems, e.g., the singlet and triplet states of spin-blockaded electrons in double quantum dots 
1 13811 . However, further work remains to be done on identifying the exact radical pair which sustains this mechanism, 
to show that it has the properties necessary to act as a compass, and to verify this with further in vitro experiments. 
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Figure 7: (Color online) (a) Tunnelling distances (Ru-Cu) and backbone structure models of several Ru-modified azurins. The azurin protein has 
a Cu active site (blue) and several His residues that can be modified to attach a Ru(bpy)2(im)(HisX)-* label (orange). This allows measurements 
of Cu(I)— >Ru(III) electron transfer at vaiious distances through the protein /?-sheet. (b) TunnelHng timetable for intraprotein electron tunnelling in 
Ru-modified proteins. In general, the distance dependence shows an exponential decay, illustrating a single-step tunnelling mechanism. The two 
solid lines depict the predictions for coupling along /3-strands (decay constant f) = l.OA"') and a-helices (fi = 1.3A"') predicted by the tunnelhng- 
pathway approach. The dashed line shows a/3 = I.IA"' dis tance decay in proteins. Distance dependence for electron tunnelling through water as 
well as vacuum are also shown. Both images are from Ref. Il39ll ) 




Figure 8: (Color online) The schematic shown above outlines the model proposed by Brookes et al l54ll to describe charge-transport-mediated 
odorant discrimination. In general, the olfactoiy receptor (G, a transmembrane protein) responds to an odorant molecule by releasing a subunit of 
a neighboring G protein, stalling a chain reaction of Ca ion influx into the cell (see the discussion of ion channels for how this process may also 
have quantum features). In the model proposed by Brookes et al [5£i, an electron source X arrives and exchanges charge with the receptor protein. 
The charge travels to a donor D in one part (helix) of the transmembrane receptor protein. The charge can then hop, or tunnel, to the acceptor A, 
in another helix of the protein, inelastically assisted by an odorant phonon emission. The charge, after arriving at D, triggers the release of the 
neighboring G protein subunit. Many precise features of this model are unknown. For example, the origin of the electron source X is unknown, but 
is conjectured to be an oxidizing agent in the cell fluid, which diffuses throughout the cell and randomly arrives at the transmembrane to initiate the 
charge-transport process. The goal of this model is to have the release of the neighboring G protein subunit depend on vibrational properties of the 
odorant molecule, and thus discriminate different molecules beyond features like size and shape. 
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4. Other possible quantum biological systems 



4.1. Tunneling in biological systems 

In many biological processes, tunneling of light particles (hydrogen or electrons) provides an alternative route to 
classical over-the-barrier reactions, resulting in new channels for chemical reactions that are forbidden in classical 
mechan ics. Tunneling in biological systems was first reported in long-range electron-transfer processes in proteins 
and hen in the transfer of a proton, hydride, or hydrogen atom in various enzymatic catalytic reactions 



1153 



Ml 



Furthermore, the simultaneous transfer of an electron and a proton fro m diff erent sites (so-called proton- 
coupled electron transfer) also play an important role in many biological functions il40ll . A significant increase in the 
enzyme catalytic efficiency or long-range electron transfer through proteins due to tunneling have been observed in 
many biological systems. This class of quantum phenomena in biology clearly fits into our description of functional 
quantum biology. 

Electron transfer between redo x cen ters separated by distances of the order of 15-30 A plays important roles in 
respiration and photosynthesis ll24[ll4lll . In the past three decades. Gray and coworkers 112411 14211 have yielded a re- 
markably detailed description of the distance- and driving-force dependencies of long-range electron tunneling rates in 
ruthenium-modified proteins. Particul arly, Ru-azurin has proven to be an excellent model system for electron tunnel- 
ing through folded proteins (Fig. 7a) 11431 . Experimentally observed electron transfer rates that exhibit exponential 
decay as distance increases and weak temperature dependence indicate a single-step tunneling mechanism (Fig. 7b). 
It is remarkable that such long-distance electron transfer in biology occurs via quantum mechanical tunneling. Tun- 
neling over such long distances would be impossible in vacuum, however, through bonded as well as non-bonded 
interactions in the protein, electron tunneling barriers can be significantly reduced, leading to significantly higher 
electron tunneling rates in proteins. 

A major open question remains in this fiel d is w hether or not the electron conducting proteins have evolved 
specific pathways for electron tunneling lll4lll24lll44ll . Theoretical analysis has demo nstrated th at an empirical model 
treating the protein as a structureless random medium explains the experimental data lll44[ll45[l . However, structure- 
based analysis reveals that specific channels through covalent bond, hydrogen bond or even van der Waa ls co ntacts 
form electron tunneling pathways that facilitate electron tunneling over a long distance 

these redox proteins specifically wired to perform efficient electron tunneling? How much do the functions of these 
biological systems depend on the quantumness of the process? Again, these are not questions that can be answered 
easily, and interpretation of current experimental results often involves ambiguities. 

Recently, Skourtis and coworkers have suggested a proposal that is likely to provide an incisive t est fo r the path- 
ways model. Quantum theory predicts that electron transfer pathways could interfere with each other lll47ll . and it has 
been suggested that electron transfer through the azurin protein depends critically by quantum interferences between 
multiple distinct pathways 114711 . Motivated by these earlier theoretical studies, Skourtis and coworkers suggested 
a molecular which-way interferometer experiment in which localized and distinguishable normal modes coupled to 
bridge ato ms can be vibrationally excited to control inelastic scattering of electrons and interferences between transfer 
pathways 1 149l 148 1 . Experiments have demonstrated that electron transfer in a a small donor- bridg e-acceptor organic 
molecule can be controlled by excitation of vibrational modes localized on the bridge moiety lIlSOll . Such experimen- 
tal probes of pathway coherence in electron-transfer proteins in the single-protein level would provide decisive proof 
for the quantumness and functionality of the electron tunneling process. 

Tunneling of a proton, hydride, or hydrogen atom also play an important role in a wide range of biological enzy- 
matic catalytic reactions i57, ,56.1. Measurements of the intrinsic kinetic isotope effects (KIE) in enzymatic reactions 
have clearly demonstrated nucl ear qu antum effects in enzymes 1 57 , 5^ . For example, in soybean lipoxygenase 1 151 1 
and methane mono-oxygenase [ISjjl, large KIEs with kn/ko close to 100 have been observed, and the sizes of the 
KIEs clearly indicate quantum-tunneling effects. In enzyme catalysis, a large portion of the quantum improvement 
over the classical catalytic rates can be attributed to the e nerg y shift due to the zero-point energy that gives a quantum 
correction to the barrier height and the H-tunneling rates pT!]. Note that semi-classical mode ls includin g environmen- 



tally coupled H-tunneling have been shown to adequately describe H-tunneling in enzymes 01531115411 . Such nuclear 



quantum effects in enzymes might represent a class of quantum phenomena in biological systems that depends only 
on trivial quantum effects, not quantum coherence. 
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4.2. Smell 

Our sense of smell allows us to discriminate between small molecules in very low concentrations via scent 
molecules interacting with receptors in the nose. At this time, the biomolecular processes of olfaction are not fully 
understood, and some evidence suggests that a mechanism based solely on the size and shape of odorant molecules 
is inadequate. For example, it has been noted that molecules with very similar shapes and sizes have a remarkably 
different scent lE6ll . Thus traditional models of a "docking"-type mechanism, where the size and shape of an odorant 
molecule actuates the receptor in some way, are thought to be insufficient. Turin [i26i1 proposed a mechanism which, in 
addition to "docking", gives a further level of selectivity (or sensitivity) by a process of inelastic electron tunnelling. 
In this case the odorant molecule both docks with a receptor and then mediates phonon-assisted inelastic tunnelling 
of an electron from a donor to an acceptor (i.e., donor and acceptor electronic states differ in energy by hu, and thus 
transport only occurs when energy is conserved by emission of an odorant phonon, the vibrational degree of freedom 
of the molecule one is "smelling", of the right energy). 

A recent model proposed by Brookes ll54ll expanded on this idea, and presented evidence that such a mechanism 
fits the observed features of smell, and is at least "physically" credible (see Fig. [S] for an overview and details of the 
mechanism). Whether such a mechanism ultimately exists in nature has yet to be determined. While not specifically 
requiring "coherence" to function, this mechanism requires inelastic phonon-assisted tunnelling of electrons, and is 
certainly more "microscopic", and sensitive, than previously thought. 



4.3. Quantum coherence in ion channels 

Ion channels, wh ich r egulate the flow of ions across the membrane of a cell, are a vital component of many 
biological processes ll55ll , e.g., neuronal communication (potassium channels), muscle contraction, etc. In many 
cases the flow of ions through the channel is controlled by a gate that can be activated by a voltage, a chemical 
signal, incident light or mechanical stress. In addition the channels often feature a filter which controls which types 
of io ns m ay pass. This selectivity filter is only a few angstroms wide, which forces the ions to pass through it one-by- 
one lfl56ll . 

A recent conjecture by Vaziri and Plenio 112911 is that the selectivity filter of ion channels, and the transport of 
ions through the channel, may exhibit quantum coherence. They are motivated by the large selectivity found in the 
bacterial KcsA potassium channel, which can select potassium over sodium with a ratio of 10^. They argue that the 
energy scale of the thermal fluctuations of the atoms which form the channel suggest two possible quantum coherent 
phenomena. Firstly, there could be a coherent diffraction of a 1-D potassium "matter wave" off an effective grating 
formed by modulations of the potential energy landscape inside the channel. Secondly, since the channel forms a 
1-D array of trapping sites, the potassium ions could take advantage of quantum tunnelling to process through the 
channel in an efficient manner They calculate a tunnelling rate which is of the same order as the decoherence time, 
and thus argue that an "environment assisted" quantum tunnelling phenomena, akin to that seen the energy transport 
in photosynthesis, is responsible for the high selectivity and throughput of these channels. However, unlike in the 
transport of energy in FMO, this process involves the movement of charges through a complex environment, with a 
correspondingly large re-organization energy, which implies the coupling to the environment may be too strong for 
coherence to thrive. 

This conjecture is intriguing, but at this time lacks experimental support (apart from the already observed re- 
markable performance of the channels). Vaziri and Plenio ll29ll propose several ways to verify the presence of this 
"functional quantum effect", focusing on collective transport resonances which, they argue, are the functionality pro- 
vided by the quantum coherent transport in this case. Furthermore, a nitrogen- vacancy (NV) probe in proximity to the 
ion channel I157il might be a promising tool to observe any possible quantum coherent effects. 



4.4. Photoreceptors 

Biological photoreceptors, i.e. light-sensing proteins, are ubiquitous in the functions of life forms on Earth. These 
proteins contain chromophores that upon excitation, undergo ultrafast chemical trans form ation on their excited-state 
potential energy surface and eventually lead to light-induced signal transduction i28lll58n . Because of the quantum 
nature of the electronically excited state and nuclear quantum coherence that often accompany excited-state molecular 
dynamics in the femtosecond time scale, quantum mechanics are required to describe these photo-activated dynamics 
in biological photoreceptors lll59ll . 
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Figur e 9: (Color online) (a) Schematic potential energy surfaces of electronic states involved in the photoisomerization of retinal (image from 
Ref. II66II ). (b) Structure model of bovine rhodopsin. The retinal motif is rendered in cyan. The quantum yield of ~ 65% is almost an order of 
magnitude higher than that observed in the absence of the protein environment. 



For example, the primary event in vision involves the photoisomerization of a retinal molecule in the trans- 
membrane protein called rhodopsin (Fig. |9), which is the most well studied and characterized photoreceptor Fol- 
lowing light absorption, this photoisomerization reaction proceeds with remarkable rate (<200 fs, one of the fastest 
chemical reactions known) and high specificity (quantum yield of ~ 0.65) IztIi . This isomerization triggers a protein 
conformational change that is subsequently amplified through protein-protein interactions with signal-transduction 
proteins, leading to the visual signal. The efficient photoisomerization process makes the retinal molecule an effective 
and reliable photo-switch. The quantum mechanical positioning of electronic states and their symmetri es enable s 
ultrafast coherent wave-packet dynamics that are responsible for the rapid photoisomerization reaction ||28[ Il58[ 1 16011 . 
In addition, the unique reactivity of rhodopsin has often been attributed to a conical intersection between its elec- 
tronic excited state and ground state lll6llll62ll . which has recently been demonstrated experimentally 115 511 . Nuclear 
quantum coherence clearly plays an important role in the remarkably rapid and efficient photo-activated one-way 
reaction dynamics in rhodopsin lll60ll . Intriguingly, coherent optical control of retinal isomerization yie lds in bac - 



teriorhodopsin, a protein closely related to vertebrate rhodopsins, has been demonstrated experimentally 11631116411 



The frontiers of the investigations of biological photoreceptors are thus to further elucidate the quantum nature that 
controls the hi gh yi eld and specificity of natural photoreceptors in order to design and direct chemical dynamics using 
coherent light 1165 j . 



5. Open issues 

5.1. Open problems in photosynthetic light harvesting 
5.1.1. Protein Environment and Spectral Density 

The protei n environments surrounding pigments in photosynthetic complexes play crucial roles in light harvest- 
ing ||I67l[168|]. On the one hand the protein scaffolding provides a rigid structure that fixes photosynthetic pigments 
in space to form an antenna network, and pigment-protein interactions also shift the transition energies of individual 
pigments to yield the energy landscape that directs energy flow towards the reaction center fi89i,i81il . These aspects of 
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proteins' roles are st atic a nd well recognized, and as a result structure-based studies have formed a large component 



in this research field 11169 1 



On the other hand, dynamical fluctuations of pigment-protein couplings are directly responsible for the energy 
transfer dynamics. However, these dynamical effects are often overlooked and treated without details, as the simple 
models of spectral densities described in Sec. 2.2. Does the popular diagonal exciton-phonon coupling model com- 
pletely capture the physical reality in photosynthetic complexes? How important (or unimportant) are off-diagonal 
exciton-phonon couplings? Are the simplified spectral density functions used in the literature adequate to describe the 
real system? Can correlated baths significantly influence exciton dynamics? And if yes, can the correlated-bath effects 
be directly observed in real systems? What are the importance of non-Markovian effects? Can we directly measure 
non-Markovianity and the memory kernel of system-bath interactions? All these questions are still unresolved and 
require a combination of experimental and theoretical studies in order to provide adequate answers. Experimentally, 
probes more specific to exciton-phonon interactions and bath dynamics must be developed. In this regard, the ri ch in - 



formation in the peak shapes and spectral dynamics of 2D electronic spectra must not be overlooked 11751 Il8l l74l 110311 



Theoretically, simulations of exciton transfer dynamics in photosynthetic systems must be expanded to inves tigate 



more general models including different types of exciton-phonon couplings and bath spectral densities 1116811 . Re- 
cently, molecular dynamics simulations combined with quantum chem istry calculations have been applied to study 
fluctuations in the FMO complex Il34ll87ll and reaction center This type of study might lead to new msights 



on the dynamical role of the protein environments and provide atomistic resolution to th e spectral densities; however, 
the interpretation of the hybrid-quantum/classical results should be treated with caution 1 171 , 170ll . 



5.1.2. Consequences of Coherence 

As we have discussed in Sec. 2, quantum coherence effects in photosynthesis have drawn intense research interest 
recently. So far, experimental and theoretical works have focused on the question of how quantum coherence effects 
contribute to improve the efficiency of energy transfer, and a measure of efficiency such as the quantum yield defined 
in Eq. (4) is often adopted as the target to benchmark the effects of quantum coherence. It is clear that site-basis 
quantum coherence, i.e. coherent delocalization of excitons, pla ys crucial roles in determining the energy landscape 



and spectral properties of photosynthetic complexes 1112111 3Lll4llI72lll73[l . However, the consequences of excitonic 



quantum coherence responsible for the reversible wave-like exciton dynamics are less clear. As we have detailed in 
Sec. 2.3, it has been speculated that excitonic coherence can improve functions such as robustness against energetic 



disorder or dynamical perturbations. In addition, Ishizaki and Fleming 113511 argued that excitonic coherence can be 
utiUzed to form energy "ratchets" for energy transfer that avoid trapping by local energy minima. However, these 
results are often model-dependent and a unifying view of excitonic coherence effects in photosynthesis seems to be 
still lacking. 

Furthermore, to establish excitonic coherence as a general and essential phenomena in photosynthetic light har- 
vesting, measurements on larger photosynthetic complexes such as LH2, LHl, and LHCII must be carried out to 
confirm the general existence of long-lasting quantum coherence. It is also necessary to consider more general initial 



conditions beyond the coherent laser-excitation condition suitable for the interpretation of laser experiments Hi 7411 . 
Does quantum coherence play a role when excitations are generated by incoherent sunlight or energy transfer from 
another antenna complex? Another open question that follows naturally after the observations of long-lasting quan- 
tum coherence in photosynthetic complexes is whether or not quantum effects (quantum coherence and entanglement) 
in photosynthetic systems can be used as a quantum resource? Here, a crucial experimental development is the ca- 
pability to prepare and measure an excitation localized on a single chlorophyll site, which would allow the quantum 
mechanical manipulation of the wavefunction of excitons. 

5.1.3. Biological Relevance 

So far excitonic quantum coherence has been observed only in iti vitro experiments, and theoretical studies have 
only established that quantum coherence could affect efficiency of energy transfer within a single complex. Howe ver, a 



connection between excitonic quantum coherence effects and biological function in vivo has yet to be established 0167 



I7IIT26I1 . In particular, the functionality of light harvesting systems depends on assemblies with a length-scale beyond 



a single complex; the whole antenna could contain tens of light-harvesting complexes. Does the single-complex 
efficiency provide enough information to determine the energy transfer efficiency of the whole antenna? Moreover, 
what is the system optimized for? It is conceivable that in the larger length scales and longer time scales that determine 
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the functionality of photosynthetic light harvesting, biological systems actually strive to achieve robustness and safety, 
instead of maximal efficiency. It is therefore necessary to investigate and model dynamics of light harvesting in larger 
assemblies including the whole antenna. Specifically, experimental probes for the combined temporal and spatial 
dynamics of exciton motions on the antenna have to be developed il76ll . In addition, evidence of long-lasting quantum 
coherence has yet to be obtained in vivo. Theoretically, models of quantum coherence effects should be considered 
for larger structures, and it is necessary to develop measures of quantum coherence that can be linked with in vivo 
observables that can be determined and tested in experiments. 



5.1.4. Design Principles 

It is surprising that after decades of research we still do not k now h ow to build artificial light-harvesting systems 
that can rival the efficiency achieved in natural photosynthesis lfl4lll77 l. We know that: (1) chromophores with strong 
absorbance must be used, (2) the antenna must contain high-density of light absorbers yet avoiding the formation 
of exciton quenching charge-transfer species, and (3) a energy gradient can help to funnel excitations towards the 
reaction center lll4ll . However, the knowledge required to implement these principles in artificial systems is lacking. 
Moreover, certainly there exist other design principles in the working of natural photosynthesis. If quantum coherence 
can help the efficiency of light harvesting, then what are the molecular architectures that can utilize this advantage to 
make a more efficient antenna? It is crucial that theoretical studies focused on how chromophores should be arranged 
are carried out and tested by experiments, maybe with the help of synthetic organic chemists, to eventually formulate 
the design principles that can help us build more efficient solar cells. 



5.2. Open problems in the radical-pair model 
5.2.1. Radical Pair Candidate Molecules 

The discovery of cryptochromes as a possible candidate to be the light-sensitive host of the magnetoreception 
radical pair mechanism gives hope to the radical-pair model. However, in-vitro experiments have yet to show that 
cryptochromes have the exquisite and anisotropic magnetic field sensitivity needed to function as a magnetic compass. 



Ideal ly cryptochromes, or other candidate molecules, need to have some of the following properties (see 1113311 and 



1 17811 for a more in-depth review of these issues): 



Spin-correlation and coherence time on the scale of geomagnetic B-field effects. The rate of coherent spin os- 
cillations induced by the geomagnetic field must exceed that of any environmental losses. A 50 jlT field induces 
oscillations on the time scale of 700 ns. This coherence, or correlation, must be protected from environmental 
losses on this time scale for significant change in the reaction rates to occur Such correlation times have been 
seen in "photolyases" radical pair systems. 

Spin-correlation time and coherence on the same time-scales as weak disruptive radio frequency fields. In 
behavioral experiments it was found that the magnetic sense was disrupted for oscillating fields with very small 
magnitudes. As with the geomagnetic sensitivity this implies the spin correlation time must survive for an 



exceptio nally long time. For example, in Ref. 04711 . Robins were found to be sensitive to fields as weak as 50 



nT. Ritz 1117811 discusses this point and argues that the effect of a 10 nT field induces a change in the reaction 



yield of 0.1% at 10 /zs, and 1% at 100 /is (a similar effect was shown in the calculation in 112311 V This implies both 
an exceptionally long spin correlation time and extreme sensitivity of the magnetoreception sense to changes in 
the reaction yield (see Fig. (|6]l). 



Optimization of the hyperfine field. Recent work by Cai et al 1113611 suggests that a highly anisotropic but weak 
nuclear hyperfine coupling can give a stronger angular dependance on the yield than a strong one. Whether this 
occurs in candidate radical pairs needs to be determined. 



5.2.2. Behavioral Tests 

There are a variety of further behavioral tests that have yet to be done il33 
stronger proof for the radical pair model: 



17811 . and which might provide even 
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Functionality window: For static fields it was found that changes in the intensity by more tha n 25% disrupted 



the magnetic compass mechanism. However, after a period of time the sense returned. Ritz 1117811 has raised 
the possibility that a behavioral test for a similar "re-adjustment over time" to the disruptive effect of oscillating 
fields will reveal if the disruption is related to the functionality window seen with static fields. 

Magnetoreception via oscillating fields: If the oscillating field amplitude is of similar intensity to the Earth's 
geomagnetic static field (in an environment where the geomagnetic field is suppressed), some of the simple 
radical pair models predict th at the oscillating field with also induce a directionally sensitive change in the 



reaction yields. Rodgers et al 1113311 argue that observation of animals responding in a directional way to this 
effect (essentially navigating by the oscillating field) would be stronger proof of the radical pair mechanism than 
the disruption effects observed so far (which in principle could be attributed to some unknown causal effect on 
the motivation of the animals). 

Apart from these directly testable issues, there remains a lot to be learned about the transduction of the radical 
pair yield to a "compass" that is biologically useful. The chain of events connecting changes in concentrations of 
molecules, into nerve impulses, and then into a signal of directionality is completely unknown. Indeed, as we have 
mentioned earlier, a large range of information is probably combined to ultimately form the final relevant directional 
cues (or instincts). This also has an impact on an understanding of why very small changes in the radical pair yields 
can apparently disrupt the magnetoreception sense. Finally, there have been pro posal s for using quantum control 
techniques, in the form of magnetic field pulses, to coherently control the spin states il79ll . If this induced a directional 
response in test animals this would be extremely difficult to explain via any alternative magnetoreception mechanism. 

Finally, a very recent experiment on pigeons has identified a neural substrate responsible for processing magnetic 
sense information 1II8OI1 . This substrate was shown to respond to the intensity, direction, and polarity of the magnetic 
field, per haps r uling out a radical pair mechanism in this case. For pigeons this is particularly interesting as other rec ent 



evidence 1112711 has contradicted the previously held notion that magnetite deposits in their beaks were responsible for 
their magnetic sense. Perhaps a third, unknown, mechanism plays a role in this particular species. Similar experiments 
on neural response in species in which the radical pair model is thought to reside (e.g., european robins) are vital. 



6. Concluding remarks 

At the atomic level, the demesne of quantum mechanics, one can always say that of course we exist in a quantum 
world (or at least a non-classical one), and that life itself is thus, in some manner, quantum. However, since the 
birth of quantum mechanics, the extent to which the features of quantum mechanics, such as quantum coherence 
or entanglement, can operate on biological time and energy scales in a "hot and wet" environment, and provide a 
biological advantage, has been seen mostly as a curiosity. 

The recent development of quantum information theory has revealed that quantum mechanics can, in some situ- 
ations, do things better than an equivalent classical theory (such as speed up computation and provide secure com- 
munication). Thus, as we have described here, the goal of "functional quantum biology" is in some sense similar; to 
enhance the efficiency of transport in a molecular chain, or t o provide an exquisite sensor of geomagnetic fields. A 
very extreme example is the idea proposed by Penrose 1II8III : that quantum computation may occur i n microtu bule 



assemblies within the neurons of the brain. This idea has attracted much attention, criticism and debate Ill82ill83ll be- 
cause it portrays the fascin ating poss ibiUty of quantum mechanical effects playing a central role in cognitive function. 



Evidence so far is lacking Ill82ill83[l . but this idea has helped stir up a broader interest in the role of quantum effects 
in biological processes, and the hunt for other candidates of functional quantum biology. 

Thus, at this time, the most promising candidates for function quantum biology remain photosynthetic units jsl 
BSSIj and magnetoreception £1, 43]. The evidence so far is in favor of both these systems not only containing 



quantum coherence in a hot and wet biological environment, but also that it is used to gain a biological advantage. 
There is a great possibility that further examples of functional quantum biology remain to be found. After all, as in a 
great number of technological breakthroughs, we usually find that nature has been there before us. 
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